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Abstract 

The MAGIC telescope is an imaging atmospheric Cherenkov telescope, designed to observe very high energy 
gamma-rays while achieving a low energy threshold. One of the key science goals is fast follow-up of the 
enigmatic and short lived gamma-ray bursts. The drive system for the telescope has to meet two basic de- 
mands: (1) During normal observations, the 72-ton telescope has to be positioned accurately, and has to track 
a given sky position with high precision at a typical rotational speed in the order of one revolution per day. 
(2) For successfully observing GRB prompt emission and afterglows, it has to be powerful enough to position 
to an arbitrary point on the sky within a few ten seconds and commence normal tracking immediately there- 
after. To meet these requirements, the implementation and realization of the drive system relies strongly on 
standard industry components to ensure robustness and reliability. In this paper, we describe the mechanical 
setup, the drive control and the calibration of the pointing, as well as present measurements of the accuracy of 
the system. We show that the drive system is mechanically able to operate the motors with an accuracy even 
better than the feedback values from the axes. In the context of future projects, envisaging telescope arrays 
comprising about 100 individual instruments, the robustness and scalability of the concept is emphasized. 

Key words: MAGIC, drive system, lACT, scalability, calibration, fast positioning 



X 
S3 



1. Introduction 

The MAGIC telescope on the Canary Island of 
La Palma, located 2200 m above sea level at 28°45' N 
and 17°54' W, is an imaging atmospheric Cherenkov 
telescope designed to achieve a low energy threshold, 
fast positioning, and high tracking accuracy [[l], 0]. 
The MAGIC design, and the currently ongoing con- 
struction of a second telescope (MAGIC II; [[si]), 
pave the way for ground-based detection of gamma- 
ray sources at cosmological distances down to less 
than 25 GeV [4] . After the discovery of the distant 
blazars lES 1218-1-304 at a redshift of 2 = 0.182 [Q 
and lES 1011-1-496 at z = 0.2\2 [6], the most re- 
cent breakthrough has been the discovery of the 
first quasar at very high energies, the flat- spectrum 
radio source 3C279 at a redshift of 2 = 0.536 
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These observational results were somewhat surpris- 
ing, since the extragalactic background radiation in 
the mid-infrared to near-infrared wavelength range 
was believed to be strong enough to inhibit prop- 
agation of gamma-rays across cosmological dis- 
tances ['8l,|9, 10]. The apparent low level of pair atten- 
uation of gamma-rays greatly improves the prospects 
of searching for very high energy gamma-rays from 
gamma-ray bursts (GRBs), cf. [11]. Their remark- 
able similarities with blazar flares, albeit at much 
shorter timescales, presumably arise from the scal- 
ing behavior of relativistic jets, the common physical 
cause of these phenomena. Since most GRBs reside 
at large redshifts, their detection at very high ener- 
gies relies on the low level of absorption [12]. More- 
over, the cosmological absorption decreases with 
photon energy, favoring MAGIC to discover GRBs 
due to its low energy threshold. 
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Due to the short life times of GRBs and the lim- 
ited field of view of imaging atmospheric Cherenkov 
telescopes, the drive system of the MAGIC telescope 
has to meet two basic demands: during normal ob- 
servations, the 72-ton telescope has to be positioned 
accurately, and has to track a given sky position, 

1. e., counteract the apparent rotation of the celestial 
sphere, with high precision at a typical rotational 
speed in the order of one revolution per day. For 
catching the ORB prompt emission and afterglows, 
it has to be powerful enough to position the tele- 
scope to an arbitrary point on the sky within a very 
short time and commence normal tracking immedi- 
ately thereafter. To keep the system simple, i.e., ro- 
bust, both requirements should be achieved without 
an indexing gear. The telescope's total weight of 
72 tons is comparatively low, reflecting the use of 
low-weight materials whenever possible. For exam- 
ple, the mount consists of a space frame of carbon- 
fiber reinforced plastic tubes, and the mirrors are 
made of polished aluminum. 

In this paper, we describe the basic properties of 
the MAGIC drive system. In section [31 the hard- 
ware components and mechanical setup of the drive 
system are outlined. The control loops and perfor- 
mance goals are described in section El while the 
implementation of the positioning and tracking al- 
gorithms and the calibration of the drive system are 
explained in section [51 The system can be scaled 
to meet the demands of other telescope designs as 
shown in section[6l Finally, in section[71and section[8l 
we draw conclusions from our experience of operat- 
ing the MAGIC telescope with this drive system for 
four years. 

2. General design considerations 

The drive system of the MAGIC telescope is 
quite similar to that of large, alt-azimuth-mounted 
optical telescopes. Nevertheless there are quite 
a few aspects that influenced the design of the 
MAGIC drive system in comparison to optical tele- 
scopes and small-diameter Imaging Atmospheric 
Cherenkov telescopes (lACT). 

Although lACTs have optical components, the 
tracking and stability requirements for lACTs are 
much less demanding than for optical telescopes. 



Like optical telescopes, lACTs track celestial ob- 
jects, but observe quite different phenomena: Opti- 
cal telescopes observe visible light, which originates 
at infinity and is parallel. Consequently, the best- 
possible optical resolution is required and in turn, 
equal tracking precision due to comparably long in- 
tegration times, i.e., seconds to hours. In contrast, 
lACTs record the Cherenkov light produced by an 
electromagnetic air-shower in the atmosphere, in- 
duced by a primary gamma-ray, i.e., from a close 
by (5 km- 20 km) and extended event with a diffuse 
transverse extension and a typical extension of a few 
hundred meters. Due to the stochastic nature of the 
shower development, the detected light will have an 
inherent limitation in explanatory power, improving 
normally with the energy, i.e., shower-particle mul- 
tiplicity. As the Cherenkov light is emitted under a 
small angle off the particle tracks, these photons do 
not even point directly to the source like in optical 
astronomy. Nevertheless, the shower points towards 
the direction of the incoming gamma-ray and thus to- 
wards its source on the sky. For this reason its origin 
can be reconstructed analyzing its image. Modern 
lACTs achieve an energy-dependent pointing resolu- 
tion for individual showers of 6' - 0.6'. These are the 
predictions from Monte Carlo simulations assuming, 
amongst other things, ideal tracking. This sets the 
limits achievable in practical cases. Consequently, 
the required tracking precision must be at least of 
the same order or even better. Although the short 
integration times, on the order of a few nanoseconds, 
would allow for an offline correction, this should be 
avoided since it may give rise to an additional sys- 
tematic error. 

To meet one of the main physics goals, the obser- 
vation of prompt and afterglow emission of GRBs, 
positioning of the telescope to their assumed sky 
position is required in a time as short as possible. 
Alerts, provided by satellites, arrive at the MAGIC 
site typically within 10 s after the outburst UM- 
Since the life times of GRBs show a bimodal dis- 



tribution nl4\ with a peak between 10 s and 100 s. To 
achieve a positioning time to any position on the sky 
within a reasonable time inside this window, i.e. less 
than a minute, a very light-weight but sturdy tele- 
scope and a fast-acting and powerful drive system is 
required. 
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Figure 1 : The top picture shows the MAGIC I telescope with the major components of the drive system. The elevation drive unit, 
from its back side, is enlarged in the bottom left picture. Visible is the actuator of the safety holding brake and its corresponding 
brake disk mounted on the motor-driven axis. The motor is attached on the opposite side. The picture on the bottom right shows 
one of the azimuth bogeys with its two railway rails. The motor is housed in the grey box on the yellow drive unit. It drives the 
tooth double-toothed wheel gearing into the chain through a gear and a clutch. 
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3. Mechanical setup and hardware components 

The implementation of the drive system relies 
strongly on standard industry components to ensure 
robustness, reliability and proper technical support. 
Its major drive components, described hereafter, are 
shown on the pictures in fig. [21 

The azimuth drive ring of 20 m diameter is made 
from a normal railway rail, which was delivered in 
pre-bent sections and welded on site. Its head is only 
about 74 mm broad and has a bent profile. The fix- 
ing onto the concrete foundation uses standard rail- 
fixing elements, and allows for movements caused by 
temperature changes. The maximum allowable devi- 
ation from the horizontal plane as well as deviation 
from flatness is ±2 mm, and from the ideal circle it 
is Ar = 8 mm. The rail support was leveled with a 
theodolite every 60 cm with an overall tolerance of 
±1.5 mm every 60 cm. In between the deviation is 
negligible. Each of the six bogeys holds two standard 
crane wheels of 60 cm diameter with a rather broad 
wheel tread of 1 10 mm. This allows for deviations in 
the 1 1.5 m-distance to the central axis due to extreme 
temperature changes, which can even be asymmet- 
ric in case of different exposure to sunlight on either 
side. For the central bearing of the azimuth axis, a 
high-quality ball bearing was installed fairly assum- 
ing that the axis is vertically stable. For the elevation 
axis, due to lower weight, a less expensive sliding 
bearing with a teflon layer was used. These slid- 
ing bearings have a slightly spherical surface to al- 
low for small misalignments during installation and 
some bending of the elevation axis stubs under load. 

The drive mechanism is based on duplex roller 
chains and sprocket wheels in a rack-and-pinion 
mounting. The chains have a breaking strength of 
19 tons and a chain-link spacing of 2.5 cm. The ini- 
tial play between the chain links and the sprocket- 
wheel teeth is about 3 mm - 5 mm, according to the 
data sheet, corresponding to much less than an arc- 
second on the telescope axes. The azimuth drive 
chain is fixed on a dedicated ring on the concrete 
foundation, but has quite some radial distance vari- 
ation of up to 5 mm. The elevation drive chain is 
mounted on a slightly oval ring below the mirror 
dish, because the ring forms an integral part of the 
camera support mast structure. 



Commercial synchronous motors (type designa- 
tion Bosch RexrothQ MHD 1 12C-058) are used to- 
gether with low-play planetary gears (type designa- 
tion alphcl^ GTS210-M02-020B09, ratio 20) linked 
to the sprocket wheels. These motors intrinsically 
allow for a positional accuracy better than one arc- 
second of the motor axis. Having a nominal power 
of 1 1 kW, they can be overpowered by up to a factor 
five for a few seconds. It should be mentioned that 
due to the installation height of more than 2200 m 
a.s.L, due to lower air pressure and consequently less 
efficient cooling, the nominal values given must be 
reduced by about 20%. Deceleration is done oper- 
ating the motors as generator which is as powerful 
as acceleration. The motors contain 70 Nm holding 
brakes which are not meant to be used as driving 
brake. The azimuth motors are mounted on small 
lever arms. In order to follow the small irregularities 
of the azimuthal drive chain, the units are forced to 
follow the drive chain, horizontally and vertically, by 
guide rolls. The elevation-drive motor is mounted on 
a nearly 1 m long lever arm to be able to compensate 
the oval shape of the chain and the fact that the cen- 
ter of the circle defined by the drive chain is shifted 
356 mm away from the axis towards the camera. The 
elevation drive is also equipped with an additional 
brake, operated only as holding brake, for safety rea- 
sons in case of extremely strong wind pressure. No 
further brake are installed on the telescope. 

The design of the drive system control, c.f. Bretz 
et al. Ill5ll . is based on digitally controlled industrial 
drive units, one for each motor. The two motors 
driving the azimuth axis are coupled to have a more 
homogeneous load transmission from the motors to 
the structure compared to a single (more powerful) 
motor. The modular design allows to increase the 
number of coupled devices dynamically if necessary, 
c.f. Bretz et al. [\^. 

At the latitude of La Palma, the azimuth track of 
stars can exceed 180° in one night. To allow for con- 
tinuous observation of a given source at night with- 
out reaching one of the end positions in azimuth, the 
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allowed range for movements in azimuth spans from 
(p = -90° to (y9 = +318°, where ip = 0° corresponds to 
geographical North, and (p = 90° to geographical 
East. To keep slewing distances as short as possi- 
ble (particularly in case of GRB alerts), the range 
for elevational movements spans from ^^ =+100° to 
6 = -10° where the change of sign implies a move- 
ment across the zenith. This so-called reverse mode 
is currently not in use, as it might result in hysteresis 
effects of the active mirror control system, still un- 
der investigation, due to shifting of weight at zenith. 
The accessible range in both directions and on both 
axes is limited by software to the mechanically ac- 
cessible range. For additional safety, hardware end 
switches are installed directly connected to the drive 
controller units, initiating a fast, controlled decelera- 
tion of the system when activated. To achieve an az- 
imuthal movement range exceeding 360°, one of the 
two azimuth end-switches needs to be deactivated at 
any time. Therefore, an additional direction switch 
is located at (/? = 164°, short-circuiting the end switch 
currently out of range. 

4. Setup of the motion control system 

The motion control system similarly uses standard 
industry components. The drive is controlled by the 
feedback of encoders measuring the angular posi- 
tions of the motors and the telescope axes. The en- 
coders on the motor axes provide information to mi- 
cro controllers dedicated for motion control, initiat- 
ing and monitoring every movement. Professional 
built-in servo loops take over the suppression of os- 
cillations. The correct pointing position of the sys- 
tem is ensured by a computer program evaluating the 
feedback from the telescope axes and initiating the 
motion executed by the micro controllers. Addition- 
ally, the motor-axis encoders are also evaluated to in- 
crease accuracy. The details of this system, as shown 
in figure [2l are discussed below. 

4. 1 . Position feedback system 

The angular telescope positions are measured by 
three shaft-encoders (type designation Hengstleilj 
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AC61/1214EQ.720LZ). These absolute multi-turn 
encoders have a resolution of 4096 (10 bit) revolu- 
tions and 16 384 (14 bit) steps per revolution, corre- 
sponding to an intrinsic angular resolution of 1.3' per 
step. One shaft encoder is located on the azimuth 
axis, while two more encoders are fixed on either 
side of the elevation axis, increasing the resolution 
and allowing for measurements of the twisting of the 
dish (fig. [3]). All shaft encoders used are watertight 
(IP 67) to withstand the extreme weather conditions 
occasionally encountered at the telescope site. The 
motor positions are read out at a frequency of 1 kHz 
from 10 bit relative rotary encoders fixed on the mo- 
tor axes. Due to the gear ratio of more than one thou- 
sand between motor and load, the 14 bit resolution of 
the shaft encoder system on the axes can be interpo- 
lated further using the position readout of the motors. 
For communication with the axis encoders, a CAN- 
bus interface with the CANopen protocol is in use 
(operated at 125 kbps). The motor encoders are di- 
rectly connected by an analog interface. 

4.2. Motor control 

The three servo motors are connected to individ- 
ual motion controller units (DKC, type designation 
Bosch Rexroth, DKC ECODRIVE03.3-200-7-FW), 
serving as intelligent frequency converters. An in- 
put value, given either analog or digital, is converted 
to a predefined output, e.g., command position, ve- 
locity or torque. All command values are processed 
through a chain of built-in controllers, cf. fig. HI re- 
sulting in a final command current applied to the mo- 
tor. This internal chain of control loops, maintaining 
the movement of the motors at a frequency of 1 kHz, 
fed back by the rotary encoders on the correspond- 
ing motor axes. Several safety limits ensure damage- 
free operation of the system even under unexpected 
operation conditions. These safety limits are, e.g., 
software end switches, torque limits, current limits 
or control-deviation limits. 

To synchronize the two azimuth motors, a master- 
slave setup is used. While the master is addressed 
by a command velocity, the slave is driven by the 
command torque output of the master. This operation 
mode ensures that both motors can apply their com- 
bined force to the telescope structure without oscil- 
lations. In principle it is possible to use a bias torque 
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Figure 2: Schematics of the MAGIC I (left) and MAGIC II {right) drive system. The sketches shows the motors, the motor-encoder 
feedback as well as the shaft-encoder feedback, and the motion-control units, which are themselves controlled by a superior control, 
receiving commands from the control PC, which closes the position-control loop. The system is described in more details in 
sectionHl 



to eliminate play. This feature was not used because 
the play is negligible anyhow. 

4.3. Motion control 

The master for each axis is controlled by preset- 
ting a rotational speed defined by ±10 V on its ana- 
log input. The input voltage is produced by a pro- 
grammable micro controller dedicated to analog mo- 
tion control, produced by Z&B0 (MACS, type des- 
ignation MACS). The feedback is realized through a 
500-step emulation of the motor's rotary encoders by 
the DKCs. Elevation and azimuth movement is reg- 
ulated by individual MACSs. The MACS controller 
itself communicates with the control software (see 
below) through a CANbus connection. 

It turned out that in particular the azimuth motor 
system seems to be limited by the large moment of 
inertia of the telescope ( ~ 4400 tm^, for com- 
parison Jci~850tm^; note that the exact numbers 
depend on the current orientation of the telescope). 
At the same time, the requirements on the elevation 
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drive are much less demanding. 

MAGIC 11 For the drive system several improve- 
ments have been provided: 

• 13 bit absolute shaft-encoders (type designation 
Heidenhairlfl ROQ 425) are installed, providing 
an additional sine-shaped ±1 Vss output within 
each step. This allows for a more accurate in- 
terpolation and hence a better resolution than 
a simple 14 bit shaft-encoder. These shaft- 
encoders are also water tight (IP 64), and they 
are read out via an EnDat 2.2 interface. 

• All encoders are directly connected to the 
DKCs, providing additional feedback from the 
telescope axes itself. The DKC can control the 
load axis additionally to the motor axis provid- 
ing a more accurate positioning, faster move- 
ment by improved oscillation suppression and 
a better motion control of the system. 
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Figure 4: The internal flow control between the individual controllers inside the drive control unit. Depending on the type of the 
command value, different controllers are active. The control loops are closed by the feedback of the rotary encoder on the motor, 
and a possible controller on the load axis, as well as the measurement of the current. 
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Figure 3: The measured difference between the two shaft- 
encoders fixed on either side of the elevation axis versus zenith 
angle. Negative zenith angles mean that the telescope has been 
flipped over the zenith position to the opposite side. The av- 
erage offset from zero corresponds to a the twist of the two 
shaft encoders with respect to each other The error bars denote 
the spread of several measurements. Under normal conditions 
the torsion between both ends of the axis is less than the shaft- 
encoder resolution. 



The analog transmission of the master's com- 
mand torque to the slave is replaced by a direct 
digital communication (EcoX) of the DKCs. 
This allows for more robust and precise slave 
control. Furthermore the motors could be cou- 
pled with relative angular synchronism allow- 
ing to suppress deformations of the structure by 
keeping the axis connecting both motors stable. 

A single professional programmable logic con- 
troller (PLC), in German: Speicherprogammier- 
bare Steuerung (SPS, type designation Rexroth 
Bosch, IndraControl SPS L 20) replaces the two 
MACSs. Connection between the SPS and the 
DKCs is now realized through a digital Profibus 
DP interface substituting the analog signals. 



• The connection from the SPS to the control PC 
is realized via Ethernet connection. Since Eth- 
ernet is more commonly in use than CANbus, 
soft- and hardware support is much easier. 

4.4. PC control 

The drive system is controlled by a standard 
PC running a Linux operating system, a custom- 
designed software based on ROOT [|l7|] and the po- 
sitional astronomy library slalib fisl] . 

Algorithms specialized for the MAGIC tracking 
system are imported from the Modular Analysis and 
Reconstruction Software package (MARS) jiol l20. 
2l|] also used in the data analysis [|22L 1231]. 



4.4.1. Positioning 

Whenever the telescope has to be positioned, the 
relative distance to the new position is calculated 
in telescope coordinates and then converted to mo- 
tor revolutions. Then, the micro controllers are in- 
structed to move the motors accordingly. Since the 
motion is controlled by the feedback of the encoders 
on the motor axes, not on the telescope axes, back- 
lash and other non-deterministic irregularities cannot 
easily be taken into account. Thus it may happen that 
the final position is still off by a few shaft-encoder 
steps, although the motor itself has reached its de- 
sired position. In this case, the procedure is repeated 
up to ten times. After ten unsuccessful iterations, the 
system would go into error state. In almost all cases 
the command position is reached after at most two or 
three iterations. 

If a slewing operation is followed by a tracking 
operation of a celestial target position, tracking is 
started immediately after the first movement with- 
out further iterations. Possible small deviations, nor- 
mally eliminated by the iteration procedure, are then 
corrected by the tracking algorithm. 
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4.4.2. Tracking 

To track a given celestial target position (RA/Dec, 
J 2000.0, FK5 [24]), astrometric and misalignment 
corrections have to be taken into account. While as- 
trometric corrections transform the celestial position 
into local coordinates as seen by an ideal telescope 
(Alt/Az), misalignment corrections convert them fur- 
ther into the coordinate system specific to the real 
telescope. In case of MAGIC, this coordinate system 
is defined by the position feedback system. 

The tracking algorithm controls the telescope by 
applying a command velocity for the revolution of 
the motors, which is re-calculated every second. It is 
calculated from the current feedback position and the 
command position required to point at the target five 
seconds ahead in time. The timescale of 5 s is a com- 
promise between optimum tracking accuracy and the 
risk of oscillations in case of a too short timescale. 

As a crosscheck, the ideal velocities for the two 
telescope axes are independently estimated using 
dedicated astrometric routines of slalib. For security 
reasons, the allowable deviation between the deter- 
mined command velocities and the estimated veloci- 
ties is limited. If an extreme deviation is encountered 
the command velocity is set to zero, i.e., the move- 
ment of the axis is stopped. 

4.5. Fast positioning 

The observation of GRBs and their afterglows in 
very-high energy gamma-rays is a key science goal 
for the MAGIC telescope. Given that alerts from 
satellite monitors provide GRB positions a few sec- 
onds after their outburst via the Gamma-ray Burst 
Coordination Network [,25,] . typical burst durations 
of 10 s to 100 s [UJO demand a fast positioning of 
the telescope. The current best value for the ac- 
celeration has been set to ll.Vmrads"^. It is con- 
strained by the maximum constant force which can 
be applied by the motors. Consequently, the maxi- 
mum allowed velocity can be derived from the dis- 
tance between the end-switch activation and the po- 
sition at which a possible damage to the telescope 
structure, e.g. ruptured cables, would happen. From 
these constraints, the maximum velocity currently in 
use, 70.4mrads^\ was determined. Note that, as 
the emergency stopping distance evolves quadrati- 
cally with the travel velocity, a possible increase of 



the maximum velocity would drastically increase the 
required braking distance. As safety procedures re- 
quire, an emergency stop is completely controlled by 
the DKCs itself with the feedback of the motor en- 
coder, ignoring all other control elements. 

Currently, automatic positioning by A(y9= 180° in 
azimuth to the target position is achieved within 45 s. 
The positioning time in elevation is not critical in 
the sense that the probability to move a longer path 
in elevation than in azimuth is negligible. Allowing 
the telescope drive to make use of the reverse mode, 
the requirement of reaching any position in the sky 
within 30 s is well met, as distances in azimuth are 
substantially shortened. The motor specifications al- 
low for a velocity more than four times higher. In 
practice, the maximum possible velocity is limited 
by the acceleration force, at slightly more than twice 
the current value. The actual limiting factor is the 
braking distance that allows a safe deceleration with- 
out risking any damage to the telescope structure. 

With the upgraded MAGIC II drive system, dur- 
ing commissioning in 2008 August, a maximum 
acceleration and deceleration of a^^ = 30 mrad s~^ 
and a2d = 90mrad/s^^ and a maximum velocity of 
fa^ = 290 mrad s"^ and = 330 mrad s^^ could be 
reached. With these values the limits of the motor 
power are exhausted. This allowed a movement of 
A(^= 1807360° in azimuth within 20s/33 s. 

4.6. Tracking precision 

The intrinsic mechanical accuracy of the tracking 
system is determined by comparing the current com- 
mand position of the axes with the feedback values 
from the corresponding shaft encoders. These feed- 
back values represent the actual position of the axes 
with highest precision whenever they change their 
feedback values. At these instances, the control de- 
viation is determined, representing the precision with 
which the telescope axes can be operated. In the case 
of an ideal mount this would define the tracking ac- 
curacy of the telescope. 

In figure [5] the control deviation measured for 
10.9 h of data taking in the night of 2007 July 22/23 
and on the evening of July 23 is shown, expressed as 
absolute deviation on the sky taking both axes into 
account. In almost all cases it is well below the reso- 
lution of the shaft encoders, and in 80% of the time it 
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does not exceed 1/8 of this value (~10"). This means 
that the accuracy of the motion control, based on the 
encoder feedback, is much better than 1' on the sky, 
which is roughly a fifth of the diameter of a pixel in 
the MAGIC camera (6', c.f. [26]). 



implementation. The astrometric corrections in use 
(fig. [6l) - performed when converting a celestial po- 
sition into the position as seen from Earth's center 
(apparent position) - take into account precession 
and nutation of the Earth and annual aberration, i.e.. 
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Figure 5: Control deviation between the expected, i.e. calcu- 
lated, position, and the feedback position of the shaft encoders 
in the moment at which one change their readout values. For 
simplicity, the control deviation is shown as absolute control 
deviation projected on the sky. The blue lines correspond to 
fractions of the shaft-encoder resolution. The peak at half of 
the shaft-encoder resolution results from cases in which one of 
the two averaged elevation encoders is off by one step. 



In the case of a real telescope ultimate limits of 
the tracking precision are given by the precision with 
which the correct command value is known. Its cali- 
bration is discussed hereafter. 

5. Calibration 

To calibrate the position command value, astro- 
metric corrections (converting the celestial target po- 
sition into the target position of an ideal telescope) 
and misalignment corrections (converting it further 
into the target position of a real telescope), have to 
be taken into account. 

5.1. Astrometric corrections 

The astrometric correction for the pointing and 
tracking algorithms is based on a library for cal- 
culations usually needed in positional astronomy, 
slalib ifisll . Key features of this library are the numer- 
ical stability of the algorithms and their well-tested 
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Figure 6: The transformation applied to a given set of catalog 
source coordinates to real-telescope coordinates. These cor- 
rections include all necessary astrometric corrections, as well 
as the pointing correction to transform from an ideal-telescope 
frame to the frame of a real telescope. 



apparent displacements caused by the finite speed 
of light combined with the motion of the observer 
around the Sun during the year. Next, the appar- 
ent position is transformed to the observer's posi- 
tion, taking into account atmospheric refraction, the 
Earth's rotation, and diurnal aberration, i.e., the mo- 
tion of the observer around the Earth's rotation axis. 
Some of these effects are so small that they are 
only relevant for nearby stars and optical astronomy. 
But as optical observations of such stars are used to 
train the misalignment correction, all these effects 
are taken into account. 



5.2. Pointing model 

Imperfections and deformations of the mechanical 
construction lead to deviations from an ideal tele- 
scope, including the non-exact alignment of axes, 
and deformations of the telescope structure. 

In the case of the MAGIC telescopes the opti- 
cal axis of the mirror is defined by an automatic 
alignment system. This active mirror control is pro- 
grammed not to change the optical axis once defined, 
but only controls the optical point spread function of 
the mirror, i.e., it does not change the center of grav- 
ity of the light distribution. This procedure is applied 
whenever the telescope is observing including any 
kind of calibration measurement for the drive system. 
The precision of the axis alignment of the mirrors is 
better than 0.2' and can therefor be neglected. 

Consequently, to assure reliable pointing and 
tracking accuracy, mainly the mechanical effects 
have to be taken into account. Therefore the tracking 
software employs an analytical pointing model based 
on the TPOINT^M telescope modeling software [|27l] . 
also used for optical telescopes. This model, called 
pointing model, parameterizes deviations from the 
ideal telescope. Calibrating the pointing model by 
mispointing measurements of bright stars, which 
allows to determine the necessary corrections, is a 
standard procedure. Once calibrated, the model is 
applied online. Since an analytical model is used, 
the source of any deviation can be identified and 
traced back to components of the telescope mount. 

Corrections are parameterized by alt-azimuthal 
terms (2^^, i.e., derived from vector transformations 
within the proper coordinate system. The following 
possible misalignments are taken into account: 

Zero point corrections (index errors) Trivial off- 
sets between the zero positions of the axes and 
the zero positions of the shaft encoders. 

Azimuth axis misalignment The misalignment of 
the azimuth axis in north-south and east-west 
direction, respectively. For MAGIC these cor- 
rections can be neglected. 

Non-perpendicularity of axes Deviations from 
right angles between any two axes in the 
system, namely (1) non-perpendicularity of 



azimuth and elevation axes and (2) non- 
perpendicularity of elevation and pointing axes. 
In the case of the MAGIC telescope these 
corrections are strongly bound to the mirror 
alignment defined by the active mirror control. 

Non-centricity of axes The once-per-revolution 
cyclic errors produced by de-centered axes. 
This correction is small, and thus difficult 
to measure, but the most stable correction 
throughout the years. 

Bending of the telescope structure 

• A possible constant offset of the mast bending. 

• A zenith angle dependent correction. It de- 
scribes the camera mast bending, which origi- 
nates by magic's single thin-mast camera sup- 
port strengthened by steel cables. 

• Elevation hysteresis: This is an offset correc- 
tion introduced depending on the direction of 
movement of the elevation axis. It is necessary 
because the sliding bearing, having a stiff con- 
nection with the encoders, has such a high static 
friction that in case of reversing the direction of 
the movement, the shaft-encoder will not indi- 
cate any movement for a small and stable rota- 
tion angle, even though the telescope is rotat- 
ing. Since this offset is stable, it can easily be 
corrected after it is fully passed. The passage of 
the hysteresis is currently corrected offline only. 

Since the primary feedback is located on the 
axis itself, corrections for irregularities of the chain 
mounting or sprocket wheels are unnecessary. An- 
other class of deformations of the telescope- support 
frame and the mirrors are non-deterministic and, 
consequently, pose an ultimate limit of the precision 
of the pointing. 

5.3. Determination 

To determine the coefficients of a pointing model, 
calibration data is recorded. It consists of mispoint- 
ing measurements depending on altitude and azimuth 
angle. Bright stars are tracked with the telescope at 
positions uniformly distributed in local coordinates. 



i.e., in altitude and azimuth angle. The real pointing 
position is derived from the position of the reflection 
of a bright star on a screen in front of the MAGIC 
camera. The center of the camera is defined by LEDs 
mounted on an ideal (±1 mm) circle around the cam- 
era center, cf. Riegel and Bretz (l^. 

Having enough of these datasets available, corre- 
lating ideal and real pointing position, a fit of the co- 
efficients of the model can be made, minimizing the 
pointing residual. 

5.3.1. Hardware and installations 

A 0.0003 lux, 1/2" high- sensitivity standard PAL 
CCD camera (type designation Watec WAT-902H) 
equipped with a zoom lens (type: Computar) is used 
for the mispointing measurements. The camera is 
read out at a rate of 25 frames per second using a 
standard frame-grabber card in a standard PC. The 
camera has been chosen providing adequate perfor- 
mance and easy readout, due to the use of standard 
components, for a very cheap price (< 500 Euro). 
The tradeoff for the high sensitivity of the camera 
is its high noise level in each single frame recorded. 
Since there are no rapidly moving objects within the 
field of view, a high picture quality can be achieved 
by averaging typically 125 frames (corresponding to 
5 s). An example is shown in figure |7l This exam- 
ple also illustrates the high sensitivity of the camera, 
since both pictures of the telescope structure have 
been taken with the residual light of less than a half- 
moon. In the background individual stars can be 
seen. Depending on the installed optics, stars up to 
12™ are visible. With our optics and a safe detection 
threshold the limiting magnitude is typically slightly 
above 9™ for direct measurements and on the order 
of 5™. . . 4™ for images of stars on the screen. 

5.3.2. Algorithms 

An example of a calibration-star measurement is 
shown in figure [8l Using the seven LEDs mounted 
on a circle around the camera center, the position of 
the camera center is determined. Only the upper half 
of the area instrumented is visible, since the lower 
half is covered by the lower lid, holding a special re- 
flecting surface in the center of the camera. The LED 
positions are evaluated by a simple cluster- finding al- 
gorithm looking at pixels more than three standard 



deviations above the noise level. The LED position 
is defined as the center of gravity of its light distri- 
bution, its search region by the surrounding black- 
colored boxes. For simplicity the noise level is de- 
termined just by calculating the mean and the root- 
mean-square within the individual search regions be- 
low a fixed threshold dominated by noise. 




Figure 8: A measurement of a star for the calibration of the 
pointing model. Visible are the seven LEDs and their deter- 
mined center of gravity, as well as the reconstructed circle on 
which the LEDs are located. The LEDs on the bottom part are 
hidden by the lower lid, holding a screen in front of the camera. 
The black regions are the search regions for the LEDs and the 
calibration star A few dead pixels in the CCD camera can also 
be recognized. 

Since three points are enough to define a circle, 
from all possible combinations of detected spots, 
the corresponding circle is calculated. In case of 
misidentified LEDs, which sometimes occur due to 
light reflections from the telescope structure, the ra- 
dius of the circle will deviate from the predefined ra- 
dius. Thus, any such misidentified circles are dis- 
carded. The radius determination can be improved 
further by applying small offsets of the non-ideal 
LED positions. The radius distribution is Gaussian 
and its resolution is cr < 1 mm (dr /r ^ 0.3%o) on the 
camera plane corresponding to ~1". 

The center of the ring is calculated as the average 
of all circle centers after quality cuts. Its resolution 
is ~2". In this setup, the large number of LEDs guar- 
antees operation even in case one LED could not be 
detected due to damage or scattered light. 



Figure 7: A single frame (left) and an average of 125 frames (right) of the same field of view taken with the high sensitivity PAL 
CCD camera used for calibration of the pointing model. The frames were taken with half moon. 



To find the spot of the reflected star itself, the same 
cluster-finder is used to determine its center of grav- 
ity. This gives reliable results even in case of satura- 
tion. Only very bright stars, brighter than 1.0™, are 
found to saturate the CCD camera asymmetrically. 

Using the position of the star, with respect to the 
camera center, the pointing position corresponding 
to the camera center is calculated. This position is 
stored together with the readout from the position 
feedback system. The difference between the tele- 
scope pointing position and the feedback position 
is described by the pointing model. Investigating 
the dependence of these differences on zenith and 
azimuth angle, the correction terms of the pointing 
model can be determined. Its coefficients are fit min- 
imizing the absolute residuals on the celestial sphere. 

5.4. Results 

Figure [9] shows the residuals, taken between 2006 
October and 2007 July, before and after application 
of the fit of the pointing model. For convenience, 
offset corrections are applied to the residuals before 
correction. Thus, the red curve is a measurement of 
the alignment quality of the structure, i.e., the point- 
ing accuracy with offset corrections only. By fitting a 
proper model, the pointing accuracy can be improved 
to a value below the intrinsic resolution of the sys- 
tem, i.e., below shaft-encoder resolution. In more 
than 83% of all cases the tracking accuracy is better 



than 1.3' and it hardly ever exceeds 2.5'. The few 
datasets exceeding 2.5' are very likely due to imper- 
fect measurement of the real pointing position of the 
telescope, i.e., the center of gravity of the star light. 

The average absolute correction applied (exclud- 
ing the index error) is on the order of 4'. Given 
the size, weight and structure of the telescope this 
proves a very good alignment and low sagging of the 
structure. The elevation hysteresis, which is intrinsic 
to the structure, the non-perpendicularity and non- 
centricity of the axes are all in the order of 3', while 
the azimuth axis misalignment is < 0.6'. These num- 
bers are well in agreement with the design tolerances 
of the telescope. 

5.4.1. Limitations 

The ultimate limit on the achievable pointing pre- 
cision are effects, which are difficult to correlate or 
measure, and non-deterministic deformations of the 
structure or mirrors. For example, the azimuth sup- 
port consists of a railway rail with some small defor- 
mations in height due to the load, resulting in a wavy 
movement difficult to parameterize. For the wheels 
on the six bogeys, simple, not precisely machined 
crane wheels were used, which may amplify hori- 
zontal deformations. Other deformations are caused 
by temperature changes and wind loads which are 
difficult to control for telescopes without dome, and 
which cannot be modeled. It should be noted that the 
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Figure 9: Distribution of absolute pointing residual on the sky 
between the measured position of calibration stars and their 
nominal position with only offset correction for both axes (red) 
and a fitted pointing model (blue) applied. Here in total 1162 
measurements where used, homogeneously distributed over the 
local sky. After application of the pointing model the residuals 
are well below the shaft-encoder resolution, i.e., the knowledge 
of the mechanical position of the axes. 



Figure 10: The distribution of mispointing measurements. The 
measurement is a direct measurement of the pointing accuracy. 
The plot shows its time-evolution. Details on the bin edges 
and the available statistics is given in the caption of table [l] 
Since the distribution is asymmetric, quantiles are shown, from 
bottom to top, at 5%, 13%, 32%, 68%, 87% and 95%. The 
dark grey region belong to the region between quantiles 32% 
and 68%. 



azimuth structure can change its diameter by up to 
3 cm due to day-night temperature differences, indi- 
cating that thermal effects have a non-negligible and 
non-deterministic influence. 

Like every two axis mount, also an alt-azimuth 
mount has a blind spot near its upward position re- 
sulting from misalignments of the axis which are im- 
possible to correct by moving one axis or the other. 
From the size of the applied correction it can be de- 
rived that the blind spot must be on the order of < 6' 
around zenith. Although the MAGIC drive system 
is powerful enough to keep on track pointing about 
6' away from zenith, for safety reasons, i.e., to avoid 
fast movment under normal observation conditions, 
the observation limit has been set to 6* < 30'. Such 
fast movements are necessary to change the azimuth 
position from moving the telescope upwards in the 
East to downwards in the South. In the case of an 
ideal telescope, pointing at zenith, even an infinitely 
fast movement would be required. 

5.4.2. Stability 

With each measurement of a calibration-star also 
the present pointing uncertainty is recorded. This al- 
lows for monitoring of the pointing quality and for 
offline correction. In figure [10] the evolution of the 



measured residuals over the years are shown. The 
continuous monitoring has been started in March 
2005 and is still ongoing. Quantiles are shown since 
the distribution can be asymmetric depending on 
how the residuals are distributed on the sky. The 
points have been grouped, where the grouping re- 
flects data taken under the same conditions (pointing 
model, mirror alignment, etc.). It should be noted, 
that the measured residuals depend on zenith and az- 
imuth angle, i.e., the distributions shown are biased 
due to inhomogeneous distributions on the sky in 
case of low statistics. Therefore the available statis- 
tics is given in tabled! 

The mirror focusing can influence the alignment 
of the optical axis of the telescope, i.e., it can mod- 
ify the pointing model. Therefore a calibration of 
the mirror refocusing can worsen the tracking accu- 
racy, later corrected by a new pointing model. Al- 
though the automatic mirror control is programmed 
such that a new calibration should not change the 
center of gravity of the light distribution, it happened 
sometimes in the past due to software errors. 

The determination of the pointing model also re- 
lies on a good statistical basis, because the measured 
residuals are of a similar magnitude as the accuracy 
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ing accuracy at least of the order of the shaft-encoder 
resolution is guaranteed. Pointing accuracy - already 
including all possible pointing corrections - is domi- 
nated by dynamic and unpredictable deformations of 
the mount, e.g., temperature expansion. 

7. Outlook 



Table 1 : Available statistics corresponding to the distributions 
shown in figure [TO] Especially in cases of low statistics the 
shown distribution can be influenced by inhomogeneous dis- 
tribution of the measurement on the local sky. The dates given 
correspond to dates for which a change in the pointing accuracy, 
as for example a change to the optical axis or the appUcation of 
a new pointing model, is known. 



of a single calibration-star measurement. The visible 
improvements and deterioration are mainly a conse- 
quence of new mirror focusing and following imple- 
mentations of new pointing models. The improve- 
ment over the past year is explained by the gain in 
statistics. 

On average the systematic pointing uncertainty 
was always better than three shaft-encoder steps (cor- 
responding to 4'), most of the time better than 2.6' 
and well below one shaft-encoder step, i.e. 1.3', in 
the past year. Except changes to the pointing model 
or the optical axis, as indicated by the bin edges, 
no degradation or change with time of the pointing 
model or a worsening of the limit given by the tele- 
scope mechanics could be found. 

6. Scalability 

With the aim to reach lower energy thresholds, the 
next generation of Cherenkov telescopes will also in- 
clude larger and heavier ones. Therefore more pow- 
erful drive systems will be needed. The scalable 
drive system of the MAGIC telescope is suited to 
meet this challenge. With its synchronous motors 
and their master-slave setup, it can easily be extended 
to larger telescopes at moderate costs, or even scaled 
down to smaller ones using less powerful compo- 
nents. Consequently, telescopes in future projects, 
with presumably different sizes, can be driven by 
similar components resulting in a major simplifica- 
tion of maintenance. With the current setup, a track- 



Currently, efforts are ongoing to implement the as- 
trometric subroutines as well as the application of 
the pointing model directly into the Programmable 
Logic Controller. A first test will be carried out 



within the DWARF project soon [12911 . The direct 
advantage is that the need for a control PC is omit- 
ted. Additionally, with a more direct communica- 
tion between the algorithms, calculating the nomi- 
nal position of the telescope mechanics, and the con- 
trol loop of the drive controller, a real time, and thus 
more precise, position control can be achieved. As 
a consequence, the position controller can directly 
be addressed, even when tracking, and the outer- 
most position control-loop is closed internally in the 
drive controller. This will ensure an even more ac- 
curate and stable motion. Interferences from exter- 
nal sources, e.g. wind gusts, could be counteracted 
at the moment of appearance by the control on very 
short timescales, on the order of milli- seconds. An 
indirect advantage is that with a proper setup of the 
control loop parameters, such a control is precise and 
flexible enough that a cross-communication between 
the master and the slaves can also be omitted. Since 
all motors act as their own master, in such a system 
a broken motor can simply be switched off or me- 
chanically decoupled without influencing the general 
functionality of the system. 

An upgrade of the MAGIC I drive system accord- 
ing to the improvements applied for MAGIC II is on- 
going. 

8. Conclusions 

The scientific requirements demand a powerful, 
yet accurate drive system for the MAGIC telescope. 
From its hardware installation and software imple- 
mentation, the installed drive system exceeds its de- 
sign specifications as given in section [2l At the same 
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time the system performs reliably and stably, show- 
ing no deterioration after five years of routine oper- 
ation. The mechanical precision of the motor move- 
ment is almost ten times better than the readout on 
the telescope axes. The tracking accuracy is domi- 
nated by random deformations and hysteresis effects 
of the mount, but still negligible compared to the 
measurement of the position of the telescope axes. 
The system features integrated tools, like an analyt- 
ical pointing model. Fast positioning for gamma- 
ray burst foUowup is achieved on average within less 
than 45 seconds, or, if movements across the zenith 
are allowed, 30 seconds. Thus, the drive system 
makes MAGIC the best suited telescope for obser- 
vations of these phenomena at very high energies. 

For the second phase of the MAGIC project and 
particularly for the second telescope, the drive sys- 
tem has been further improved. By design, the drive 
system is easily scalable from its current dimensions 
to larger and heavier telescope installations as re- 
quired for future projects. The improved stability is 
also expected to meet the stability requirements, nec- 
essary when operating a larger number of telescopes. 
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